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)£ ' Abstract 

ON ' 

r-j ■ By using a generalization of the optical tomography technique we describe 

Qh! the dynamics of a quantum system in terms of equations for a purely classical 

probability distribution which contains complete information about the system. 

a 

1 Introduction 
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Due to the Heisemberg |T| and Schrodinger- Robertson 0, uncertainty relation for 

the position and momentum in quantum systems, does not exist joint distribution 

function in the phase space. Nevertheless, a permanent wish to understand quantum 

mechanics in terms of classical probabilities leads to introduce the so called quasi- 

probability distributions, such as Wigner function Q, Husimi Q-function § and 

Glauber-Sudarshan P-function f|, @. Later on a set of s-ordered quasi-distributions 

p| unified these quasi-probabilities into one-parametric family. Even in the early 

days of quantum mechanics Madelung || observed that the modulus and the phase 
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of wave function obey the hydrodinamical classical equations, and along this line the 
stochastic quantization scheme has been suggested by Nelson to link the classical 
stochastic mechanics formalism with the quantum mechanical basic entities, such as 

was 



wave function and propagator. In some sense, also the hidden variables [11 



proposed to relate the quantum processes to the classical ones. Nevertheless, up to 
date there not exsist a formalism which consistently connects the "two worlds". 

The discussed quasi-probabilities illuminated the similarities and the differences 
between classical and quantum considerations, and they are widely used as instrument 



for calculations in quantum theory []T2| , ||13| . However, they cannot play the role of 
classical distributions, since for example, the Wigner function and the P-function may 
have negative values. Althought the Q-function is always positive and normalized, it 
does not describe measurable distributions of concrete physical variables. 

Using the formalism of Ref . |§ , Vogel and Risken 0] found an integral relation 
between the Wigner function and the marginal distribution for the measurable ho- 
modyne output variable which represents a rotated quadrature. This result gives the 
possibility of measuring the quantum state, and it is referred as optical homodyne 
tomography fl5| . 

In Ref. |L6| a symplectic tomography procedure was suggested to obtain the 
Wigner function by measuring the marginal distribution for a shifted and squeezed 
quadrature, which depends on extra parameters. In Ref. []T7| the formalism of Ref. 
|Tj| was formulated in invariant form, relating the homodyne output distribution 



directly to the density operator. In Ref. |18| the symplectic tomography formalism 
was also formulated in this invariant form and it was extended to the multimode 
case. Thus, due to the introduction of quantum tomography procedure the real 
positive marginal distribution for measurable observables, such as rotated shifted 
and squeezed quadratures, turned out to determine completely the quantum states. 
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The aim of the present work is to formulate the standard quantum dynamics in 
terms of the classical marginal distribution of the measurable shifted and squeezed 
quadrature components, used in the symplectic tomography scheme. Thus we obtain 
an alternative formulation of the quantum system evolution in terms of evolution of 
real and positive distribution function for measurable physical observables. We will 
show the connection of such "classical" probability evolution with the evolution of 
the above discussed quasi-prbability distributions. 

Examples relative to states of harmonic oscillator and free motion will be consid- 
ered in the frame of the given formulation of quantum mechanics. 

2 Density operator and distribution for shifted and 
squeezed quadrature 

In Ref. [[L£| it was shown that, for the generic linear combination of quadratures, 
which is a measurable observable {% = 1) 

X = nq + up + 5, (1) 

where q and p are the position and momentum respectively, the marginal distribution 
w(X, (/,, u, S) (normalized with respect to the X variable), depending upon three extra 
real parameters //, u, 5, is related to the state of the quantum system, expressed in 
terms of its Wigner function W(q,p), as follows 

w(X,p,,,5) = J e-^ x ^-^W(q,p)^^. (2) 

This formula can be inverted and the Wigner function of the state can be expressed 



in terms of the marginal distribution [16 



W(q,p) = (2irys 2 e™*w F (X,sq,sp,s), (3) 
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where Wp(X, a, b, s) is the Fourier component of the marginal distribution (fj) taken 
with respect to the parameters \i, v, 5, i.e. 

w F (X, a, b, s) = I w(X, /I, v, 5)e-^ a+vb+&s) dpdvd5. (4) 

(27r) d J 

Hence, it was shown that the quantum state could be described by the positive 
classical marginal distribution for the squeezed, rotated and shifted quadrature. In 
the case of only rotated quadrature, p = cos0, v = sin0 and 5 = 0, the usual 



optical tomography formula of Ref. |T3, gives the same possibility through the 



Radon transform instead of the Fourier transform. This is, in fact, a partial case of the 
symplectic transformation of quadrature since the rotation group is a subgroup of the 
symplectic group ISp(2, R) whose parameters are used to describe the transformation 
©■ 

In Ref. [18| an invariant form connecting directly the marginal distribution 
w(X, n, v, S) and the density operator was found 

p = J dpdvdb w(X, p,v,5)K^ U) s, (5) 

where the kernel operator has the form 

The formulae (j^) and (|5p of symplectic tomography show that there exist an invertible 
map between the quantum states described by the set of nonnegative and normalized 
hermitian density operators p and the set of positive, normalized marginal distribu- 
tions ("classical" ones) for the measurable shifted and squeezed quadratures. So, the 
information contained in the marginal distribution is the same which is contained in 
the density operator; and due to this, one could represent the quantum dynamics in 
terms of evolution of the marginal probability. 
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3 Quantum evolution as classical process 

We now derive the evolution equation for the marginal distribution function w using 
the invariant form of the connection between the marginal distribution and the density 
operator given by the formula @. Then from the equation of motion for the density 
operator 

d t p = -t[H,p] (7) 
we obtain the evolution equation for the marginal distribution in the form 

J dpdudd \w(X, fi, v y 5, t)K^ & + w(X, /i, u, 5, t)I^s\ = (8) 

in which the known Hamiltonian determines the kernel I^vs through the commutator 

i^ s = i[H,k^ s \. (9) 

The obtained integral-operator equation for simple cases can be reduced to the partial 
differential equation. To do this we represent the kernel operator 1^,8 in normal order 
form (i.e. all the momentum operators on the left side and the position ones on the 
right side) containing the operator K^s as follow 

: 4„, a := Tl{p) : K^ s : V(q) (10) 

where TZ(p) and V(q) are, finite or infinite operator polynomials (depending also 
on the parameters \x and v) determined by the Hamiltonian. Then calculating the 
matrix elements of the operator equation (H) between the states (p\ and \q) and using 
the completness property of the Fourier exponents we arrive at the following partial 
differential equation for the marginal distribution function 

d t w + U(p, q)w = (11) 
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where the polynomial Tl(p, q) is the product of the polynomials TZ(p) and V(q) rep- 
resented in the form 



U(p,q) = K(jp)V{q) =EEA%™(ft^) (12) 

n m 

in which the c-number variables p and q should be replaced by the operators 
,_ ( 1 d .p,d_\ ~_(_}_^_ 

P ~ {d/d5dv +l 2d5)' q ~\d/d8dv l 2d8) ] { } 

where derivative in the denominator is understood as integral operator. One should 
point out that the operators p and q in Eq. ([11]) act on the product of coefficients 
c n m (yU, z/) and the marginal distribution corresponding to the order shown by Eqs. 
( |TI| ) and flT2|). Let us consider the important example of the particle motion in a 
potential with the Hamiltonian 

H = P - + V{q)- (14) 

then the described procedure of calculating the normal order kernel ( |i~0|) gives the 
following form of the quantum dynamics in terms of a Fokker-Planck-like equation 
for the marginal distribution 



d 

w — u-^-w — I 
ov 



'( 1 d v d\ / 1 d ud 

V 7777777^ +«XTT7 ~ V[ 77777777" ~ %: 



w = (15) 



d/dSdfx 2 85) \d/85d^ 2d5 / 

which in general case is an integro-differential equation. For the free motion, V = 0, 
this evolution equation becomes the first order partial differential equation 

d 

w — W77-W = 0. (16) 
ov 

For the harmonic oscillator, V(q) = q 2 /2, the quantum dynamic equation has the 
form 

3d 

w — u— w + v-— w = 0. (17) 

OV Ojl 



Thus given a Hamiltonian of the form ( JTjj ) we can study the quantum evolution of 



the system writing down a Fokker-Planck-like equation for the marginal distribution. 
Solving this one for a given initial positive and normalized marginal distribution we 
can obtain the quantum density operator p(t) according to Eq. (H). Conceptually it 
means that we can discuss the system quantum evolution considering classical real 
positive and normalized distributions for the measurable variable X which is shifted 
and squeezed quadrature. The distribution function which depends on extra param- 
eters obeys a classical equation which preserves the normalization condition of the 
distribution. In this sense we always can reduce the quantum behaviour of the sys- 
tem to the classical behaviour of the marginal distribution of the shifted and squeezed 
quadrature. Of course, this statement respects the uncertainty relation because the 
measurable marginal distribution is the distribution for one observable. That is the 
essential difference (dispite of some similarity) of the introduced marginal distribution 
from the discussed quasi-distributions, including the real positive Q-function, which 
depend on the two variables of the phase space and are normalized with respect to 
these variables. We would point out that we do not derive quantum mechanics from 
classical stochastic mechanics, i.e. we do not quantize any classical stochastic pro- 
cess, our result is to present the quantum dynamics equations as classical ones, and in 
doing this we need not only classical Hamiltonian but also its quantum counterpart. 

4 Examples 

Below we consider simple examples of the marginal distribution evolution for states 
of free motion and harmonic oscillator. First of all we take into account the free 
motion for which the Eq. flTBD has a gaussian solution of the form 
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where the dispersion of the observable X depends on time and parameters as follow 

1 



= |[/i 2 + ^ 2 (l + t 2 ) + 2^t]. 



(19) 



The initial condition corresponds to the marginal distribution of the ground state of 
an artificial harmonic oscillator calculated from the respective Wigner function [|nj . 
If we consider the first excited state of the harmonic oscillator, we know the 



Wigner function |19 



Wi(g,p) = -2(1 -2q 2 - 2p 2 ) exp[-g 2 - p 2 



(20) 



It results time independent due to the stationarity of the state, but for small q and 
p it becomes negative while the solution of Eq. ([H]) 



Wl (X,n, v ,S,t) = + v 2 ]-i{X - 5) 2 exp {-^r-^r 



(21) 



is itself time independent, but everywhere positive. 

Indeed, a time evolution is present explicitly in the coherent state, whose Wigner 
function is given by 

W c (q,p) = 2exp{-g 2 - q 2 - p 2 - pi + 2(qq + pp ) cost - (pq - qp ) sint} (22) 

where go an d Po are the initial values of position and momentum. For the same state 
the marginal distribution shows a more complicate evolution 

1 



w c (X,fi,u, 8,t) = — =[// +v 



21-7 
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(23) 



(X-5) 2 (X-S), 
x exp <; -g - Po — h 2 (p cos t - g sin t) 



v 



x exp 



/i 2 + v 1 



— {X — 5) + g (/(xsint + z/cost) + p (u sint — //cost) 
v 



It is also interesting to consider the comparison between Wigner function and marginal 
probability for non-classical states of the harmonic oscillator, such as female cat state 
defined as [^0 

\ a _) = M(| a > - | - a)), a = 2" 1 / 2 (g + ip Q ) (24) 
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with 

_( exp[(g 2 +pg)/2] \* 
~ l4sinh[( % 2 +^)/2]/ 

and for which the Wigner function assumes the following form 



(25) 



W-(q,p) = 2iV 2 e q2 p2 { e g ° p o cosh[2(gg + PPo) cost + 2(qp - pq ) sint] 

- cos[2(gpo — PQo) cost — 2(qq +ppo) sint]}. (26) 

The corresponding marginal distribution is 

w-(X,n,v,5,t) = Nl[wA(X,fx,u,8,t) - w B (X,n,v,5,t) 

- w B (X,/i,u,8,t) -w A (-X,/i,u,-8,t)] (27) 



with 



w A (X, /I, u, 8, t) 



+ z, 2 



(28) 



(X-8) 2 JX-8), . ' 

x exp <; -q -p — h 2 (p cos t - q sin t) 



v 



x exp | 



— (X — <5) + 5o(/isint + z/cost) + p (* / sint — /icost) 



and 



% (X, /l, I/, 5, t) = ^[/i 2 + z/ 2 ] 3 

'7T 



(29) 



(X-5) 2 n .(X-<5), . ' 

x exp<;-g -p — 22 (q cos t + p sin i) 



/ -1 



— i— (X — <5) + q (ficost — usint) + p (fisint + ucost) 



v 



The presented examples show that for the evolution of the state of a quantum system, 
one could always associate the evolution of the probability density for the random 
classical variable X which obeys "classical" Fokker-Planck-like equation, and this 
probability density contains the same information (about quantum system) which 
is contained in any quasi-distribution function. But the probability density has the 
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advantage to behave completly as the usual classical one. The physical meaning of 
the "classical" random variable X is transparent, it is considered as the position in 
an ensamble of shifted, rotated and scaled rest frames in the classical phase space of 
the system under study We could remark that for non normalized quantum states, 
like the states with fixed momentum (De Broglie wave) or with fixed position, the 
introduced map in Eq. (|5]) may be preserved. In this context the plane wave states of 
free motion have the marginal distribution corresponding to the classical white noise. 

5 Conclusions 

We have shown that it is possible to bring the quantum dynamics back to classical 
description in terms of a probability distribution containing (over) complete informa- 
tion. The time evolution of a measurable probability for the discussed observables 
could be useful both for the prediction of the experimental outcomes at a given time 
and, as mentioned above, to achieve the quantum state of the system at any time. 
Furthermore the symplectic transformation of Eq. (|l|) could be represented as a com- 
position of shift, rotation and squeezing. So, the measurement of a shifted variable 
means the measure of the coordinate in a frame in which the zero is shifted. This 
could be implemented for example by measuring the oscillator coordinate using an 
infinite ensamble of frames which are shifted with respect to the initial one (related 
method was discussed also in Ref. pl|). Furthermore if one considers the variable q 
as the photon quadrature, which corresponds to the amplitude of the electric vector 
vibrations, a rotation means a homodyne measurement, while the squeezing means 
measurement after amplification or attenuation. So, we would emphasize that our 
procedure allows to transform the problem of quantum measurements (at least for 
some observables) into a problem of classical measurements with an ensamble of 
shifted, rotated and scaled reference frames in the (classical) phase space. 
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We also want to remark that in some situations the measurements of instataneous 
values of the marginal distribution for different values of the parameters is replaced by 
measuring the distribution for these parameters which evolve in time. Such measure- 
ments may be consistent with the system evolution if the parameters time variation 
is much faster than the natural evolution of the system itself. In this case the state 
of the system does not change during the measurement process and one obtains the 
instant value of the marginal distribution and of the Wigner function. 

Finally we belive that our "classical" approach could be a powerful tool to investi- 
gate complex quantum system as for example chaotic systems in which the quantum 
chaos could be considered in a frame of equations for a real and positive distribution 
function. 
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